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Abstract 
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I. INTRODUCTION 

Recently, in a very interesting paper, Zheng et al., 1] have studied the collective di- 
rected transport of symmetrically coupled lattices in symmetric periodic potentials. They 
show that under the action of an external wave, that breaks the spatiotemporal symmetry 
and introduces inhomogeneities of the lattice, a net unidirectional current can be observed. 
Apparently the current originates from the collaboration of the lattice and the wave (am- 
plitude, frequency, and phase shifts). The study of directed transport of particles is very 
important in the physics of molecular motors 3, [1 vortex dynamics in superconductors 
J, 5j|, Josephson junction lattices nanotechnology 0,0] and many other systems. 

Many studies have been dedicated to directed transport in spatiotemporal systems 0, 
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11 



15 



stochastic forces 113 



fisll . These include systems with deterministic ac drives and 



Very important are models composed of a symmetrically coupled lattice in a symmetric 
potential field, which is driven by an external wave [jj. 
The equation studied by Zheng et al., jl] is the following: 



8j = -d sin 6j + k (9 j+1 - 29j + 6^) + A cos (cot + j<f>) . 



(1) 



This is a Frenkel-Kontorova model driven by the wave A cos (cot + 

The average current introduced by them J = Uirt^nt Ylf=i Jo (t)dt was shown to be 
nonzero for certain values of the parameters of the system. 

We wish to re-state this problem in terms of the directed motion of solitons. Solitons 
are considered as particles in many systems. Moreover, in many systems, the solitons are 



as particles 
QHQ, 



22 



231 ] . So it is very important to study the directed transport 



charge carriers 
of these objects. 

In the present paper we investigate the generalized Klein-Gordon equation: 



4>tt + l<Pt - <Pxx - G(4>) = F(x - wt), 



(2) 



where G(4>) = —dU(<p)/d<p, U(4>) is a symmetric potential with, at least, two minima at 
<p = 0i, (j) = 03, and a maximum at = 02 (0i < 02 < 03); F(x — wt) describes an external 
wave moving with velocity w. Both the sine-Gordon and 4 equations are particular cases 
of equation (J2J). The sine-Gordon equation is a continuous relative of Eq. (JTJ). 
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These equations possess important applications in condensed matter physics. For in- 
stance, in solid sate physics, they describe domain walls in ferromagnets and ferroelectric 



materials, dislocations in crystals, charge-density waves, interphase boundaries in meta 



24, 
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al- 

m. 



loys, fluxons in long Josephson junctions and Josephson transmission lines, etc 

Our initial condition will be a kink-soliton whose center-of-mass is situated at rest in 
the point x = 0. We will show that the transport of the soliton by the wave depends on the 
shape of the wave and the parameters of the system. 

We will also investigate the stabilization and transport of bubbles using external pertur- 
bations in 3.D-systems as the following 

( f )t -V 2 (f ) -G(<j ) ) = F(x } y,z,t). (3) 

We will also present the results of real experiments with laser-induced bubbles in liquids. 

The bubbles can be trapped, stabilized and transported by a laser beam of relatively low 
power. We will show that the directed transport of solitons and bubbles are related. For 
both phenomena some conditions for the external perturbations must be satisfied. 



II. STABILITY AND DYNAMICS OF SOLITONS 



Before discussing equation © we would like to present some results about the dynamics 

nnn 

of solitons in the presence of time-independent inhomogeneous perturbations (23, 12a, 29, 
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321 ] as in the following equation: 



<l>tt + l<t>t-(l>x X -G{4>) = F{x). 



(4) 



Suppose we are interested in the stability of a soliton situated in equilibrium positions 
created by the inhomogeneous force F(x). We construct a solution <j>k(x) with the topological 
and general properties of a kink. Then, we solve an inverse problem such that F(x) possesses 



the properties of the physical system we are studying 
investigate the stability of the solution. 

For this we should solve the spectral problem 
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0,0,0, alQ- 



Now we can 



where L = — d rr — 



ggM) 



4>=<t>k{x) 



Lf(x) = Tf(x), 
, t) = cj> k {x) + f(x)e xt , T = -(A 2 + 7 A). 



(5) 



The results obtained with this function can be generalize d q ualitatively to other systems 
topologically equivalent to the exactly solvable systems 
Let us see the following example: 



*>u + K<pt ~ <Pxx + sin <p = Fi(x) 



(6) 



where = 2{B 2 - l)sm\i{Bx)/ cos\i 2 {Bx). 

The exact stationary solution of Eq. © is 4>k{x) = 4 arctan[exp(i?x)]. This solution 
represents a kink-soliton equilibrated on the position x = 0. The stability problem (0) can 
be solved exactly. The eigenvalues of the discrete spectrum are given by the formula 



T„ = B 2 (A + 2An 



n 



(7) 



where A(A + 1) = 2/B 2 . The integer part of A, i.e., [A], yields the number of eigenvalues in 
the discrete spectrum, which correspond to the soliton modes (this includes the translational 
mode r , and the internal or shape modes T n (n > 0)). 

Everything that happens with this solution near point x = can be obtained from Eq. 0. 

Now let us return to the general case (J3J) and a review of the obtained results about the 
dynamics of solitons in Eq. (@J) . 

Wherever F(x) is positive, the kink-soliton will be accelerated to the "left". Wherever 
F(x) is negative, the kink-soliton will be accelerated to the "right". The zeroes of F(x) 



are candidates for equilibrium positions for the kink-soliton [27|, |28|, |29|, |30|, |31|, |32j. If 
F(x) possesses only one zero Xq (F(xq) = 0), then it is a stable position for the soliton if 
fglj.* > 0. Otherwise, the position is an unstable equilibrium. The opposite is true for the 
antikink-soliton. A soliton can be trapped by a F(x) with a zero that represents a stable 
position, and it can move away from an unstable position. 

Now we will discuss what happens in the Eq. © with different waveforms. 



III. TRANSPORT OF SOLITONS USING EXTERNAL PERTURBATIONS 

Fig. ^ shows different waveforms for F(x — vt). In case a) the kink will move trapped 
inside the potential well created by F(x — vt). The effective potential V(x) of this force F(x) 
possesses the property that V(x) — > oo when x — > ±oo. Thus, the kink will move inside a 
potential equivalent to a harmonic oscillator. So a wave like this will always transport the 
soliton to the "right" . 
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FIG. 1: Different waveforms for the transport of solitons in Eq. (J2J). a) The soliton feels a moving 
deep stable potential well, b) The soliton feels an unstable equilibrium position, c) The soliton 
feels a potential well with finite walls, d) The soliton feels a potential with a finite maximum. 

In the case b), the wave will "push" the soliton to the "right" as long as the center-of- 
mass of the soliton is situated righter than the zero of F(x). In fact, V(x) has an absolute 
maximum. However, there is a limit for the wave velocity. If the wave velocity exceeds 
this limit, the center-of-mass of the soliton can reach the maximum of the potential V(x), 
and the center-of-mass of the soliton will be behind the maximum of the potential. In this 
situation the wave can no longer push the soliton to the "right" . 

We will here introduce the concept of "controlled transport" when the soliton can be 
transported by the external wave and the soliton velocity can be approximately equal to 
the wave velocity. In the case a) of Fig. ^ we have a controlled transport. In the case 
b), the soliton is pushed by the wave but the soliton will move ahead of the wave in an 
"uncontrolled" fashion. 

The waveform shown in Fig. He represents a force that creates a stable equilibrium for the 
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FIG. 2: Effective potential "created" by the forces shown in Fig. ^ (c) and (d) respectively. 

soliton. So, as in case a), this wave can carry the soliton. However, the effective potential 
corresponding to this force is finite for x — > ±00 (See Fig. |2]a). 
So here there is also a limit for the wave velocity. 

The wave represented in Fig. dd is an unstable equilibrium for the soliton as in case b). 
This wave can push the soliton to the "right" provided its velocity is smaller that certain 
finite value. 

Fig. E] shows some different localized wave forms. 

We would think that these waves could also be used to transport a soliton in a controlled 
fashion because there is a stable equilibrium involved (provided the wave velocity is smaller 
than certain limit). However, here there are some additional limitations in the shape of the 
wave. 

Investigating different functions in (j3J) and the stability problem (|5|L we have found that 
if F(x) possesses several zeroes and the distance between them is smaller than the soliton 
width, surprising phenomena can occur. For instance, in the presence of forces as that shown 
in Fig. El a, if the zeroes of F(x) are too close, the center zero is not a stable point anymore. 
In practice, the soliton does not "feel" a minimum there. It feels an unstable equilibrium 
point. In this case, this wave would not be able to transport the soliton in a controlled 
fashion, but it can push the soliton to the right under certain conditions. 

On the other hand, the force shown in Fig. El c can carry the soliton if the zeroes are 
sufficiently separated. Otherwise, the soliton would not "feel" the left maximum of the 
potential V(x) (Fig. El d). So, using this wave the soliton cannot move trapped in the 
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FIG. 4: Perturbation equivalent to a wall. 

minimum. 

Even localized forces without a zero like that represented in Fig.|3]can "push" the soliton 
to the "right" . 
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FIG. 5: Sinusoidal wave. 



Note that, considering the effective potential V(x), this wave "acts" on the soliton as a 
moving "wall" . But once and again, there is a limit for the velocity of the wave. 

Later we will present several examples of systems of type (J2J) for which we can calculate 
all these limit velocities. 

Before, we should discuss the most ubiquitous kind of waves: the sinusoidal wave (See 



The previous discussion of waves with the shape shown in Fig. |3] leads us to think that 
the parameters A, k and w should satisfy some conditions for a soliton to be transported by 
wave ((H). In fact, if the velocity w is larger than some critical value, the soliton will be left 
behind. If k is too large (and A is small), the distance between the zeroes of F(x) will be 
very small and the potential wells, where the soliton can be trapped, would be so shallow 
that, already with certain small velocity w, the soliton will be left behind. Thus there are 
waves of type (JH) that cannot carry the soliton. 

Let us present here some concrete examples. 

Consider the following perturbed 4 equation: 



where F(x—wt) = Atanh[B(x—wt)], A > 0, B > 0. As can be expected from our discussion, 
this external wave is able to transport a soliton to the right in a controlled fashion. This 
can be observed in Fig. El 



Fig. ED: 



F(x — wt) = A sm[k(x — wt)]. 



(8) 




(9) 
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FIG. 6: The soliton is transported by the moving perturbation F{x — wt) = Atanh[B(x — wt)} 
(See Eq. @). 

Later we will see that F(x — wt) can be in practice not only a wave but a controlling 
external perturbation. If this external perturbation can be manipulated at wish in such a 
way that the velocity w(t) can be any function of time that we desire, then we can move 
the soliton and place it wherever we wish. 

Suppose now that in Eq. (JHJ) F(x — wt) = — Ata.nh[B(x — wt)}, A > 0, B > 0. In this 
case, the soliton is pushed to the "right" while it is in the "zone" where F(x — wt) < (i.e., 
while the center-of-mass of the soliton is situated righter than the zero of F(x — wt).) Let 
us define 0i < 02 < 03 as the roots of the algebraic equation — 3 = 2A. The maximum 
velocity v m that the solution can reach in its motion to the right |28| can be calculated using 
the equation 

"" (10) 



^< 2 

If the external wave velocity is such that w > v m , then the soliton eventually will start 
to move backward. 

Now we will analyze the following wave F(y), where y = x — wt: 



F(y) = Fi(y), for y > y* 
F(y) = c, for y < y* 



(11) 
(12) 
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where F x {y) = \At&vh(By)[(A 2 - 1) + (AB 2 - A 2 )sech 2 (By)), y* is the point where F x (y) 
has a local maximum; and c = Fi(y*). 

This problem is equivalent to that shown in Fig. 01c and Eld provided A 2 > 1, AB 2 < 1. 
The effective potential possesses a local minimum and a local maximum at the left side of 
the minimum. 

Using the results of papers we can solve exactly the dynamical problem of the 

soliton in this potential. 

Our analysis shows that if 2B 2 (3A 2 — 1) < 1, then the soliton feels the barrier in the 
point y = 0. So, if the wave velocity w is small, the soliton can be carried to the right. 
On the other hand, if 2B 2 {3A 2 — 1) > 1 (although the waveform is still as that shown in 
Fig-Elc such that F(y) has two zeroes), the soliton will move to the left, crossing the barrier 
even if its center-of-mass is placed in the minimum of the potential and its initial velocity 
is zero. The soliton does not feel the barrier. Even with an infinitesimal wave velocity, the 
soliton will not move to the right dragged by the wave. We have checked this phenomenon 
numerically. 

The effect of a wave shaped as that in Fig. 0] can be seen in the following example with 
exact solution: 

11 A 

kt + i<Pt - <p xx - -<P + t;<P 3 = — — ■ (13) 



- - cosh 2 



1 / x—wt 



2 \ VT- 



•w- 



When A = —\ , w 2 , the following travelling soliton is an exact solution 



tanh 



1 / x—wt 



2 VT- 

- } , i . Note that a wave shaped like that in Fig. 0] can make the soliton move 
to the right. However, for a fixed A in Eq. (p~3|) . the wave velocity should not exceed a limit. 
Otherwise, the soliton can be left behind. 

In the case of a sinusoidal wave acting in the framework of the (j) 4 equation: 

<Ptt + 70< - 4>xx - + ^<p 3 = A sin [k (x - wt)] , (14) 

we can write the following conditions for the soliton to be dragged by the wave to the right: 
t > 4, w < v m , where v m is defined by the equation 1 . Vm = and <b 2 is the "middle" 
root to the algebraic equation <fi — <p z = 2 A as in Eq. ((7j). 

We have experimented with other waveforms that illustrate the different behaviors dis- 
cussed above. 
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FIG. 7: A kink-soliton can be pushed by unstable equilibrium. F(x,t) = atanh[B(x — wt)], 
(j)(x,0) = 2Btanh[B(x-x Q )], &(s,0) = 0, B = 0.5, w = 0.05, a = -0.05, x = 1.0, 7 = 0.5 




FIG. 8: The soliton is transported by the perturbation: F(x,t) = asinh[l?(x — wt)]/ cosh 3 [i?(x — 
wt)], </>(x,0) = 2Btanh(Bx), t (x,O) = 0, B = 0.5, w = 0.1, a = 0.2, 7 = 0.5 

The force F(x,t) = atanh[S(x — wt)], with a < has the properties of the function 
shown in Fig. ^p. 

Fig. [7| shows that the soliton can be pushed by this external perturbation, provided the 
velocity w is not too large. 

In Fig. IS] we present the behavior of the soliton under the action of the external pertur- 
bation Fix, t) = as ™a!f < f~ wt )} , which is of the kind shown in Fig. [Tb. 

v 1 ' cosh [B(x— wt)\ 

When w = 0.1, the soliton is transported in a controlled way by the external wave. 
However, already when w = 0.5, the soliton is abandoned behind. This can be observed 
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FIG. 9: The same experiment as in Fig. |HJ However, the wave velocity is larger, and the soliton is 
not transported by the wave. (F(x,t) = a sinh[£>(x— wt)]/ cosh i [B(x— wt)], 4>(x,0) = 2B tanh(i?x), 
<j) t (x,0) = 0, B = 0.5, w = 0.5, a = 0.2, 7 = 2.0) 
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FIG. 10: A periodic wave can carry the soliton. ( F(x,t) = Asiii[kx — wt], (f)(x,0) = 2B tanh(i?x), 
(f>t(x,0) = 0, B = 0.5, w = 0.01, A = 0.05, k = 0.2 7 = 0.5) 

in Fig. El 

A periodic wave like F(x,t) = Asin[kx — wt] is able to carry the soliton as is shown 
in Fig. [HH where w = 0.01. The values of the other parameter can be found in the figure 
caption of Fig. [TUl 

If the wave velocity is too large (say w = 0.9), then the soliton cannot be carried by the 
external wave (See Fig. [TTJ). 
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FIG. 11: The periodic wave must satisfy some combinations in order to transport the soliton. 
In this case, for a larger wave velocity, the soliton is left behind. (F(x,t) = Asin[nx — wt], 
0(x,O) = 2Btanh(Bx), <f> t (x,0) = 0, B = 0.5, w = 0.9, A = 0.05, k = 0.2 7 = 0.5) 



IV. BUBBLES IN SOLITON EQUATIONS 

34, 35, 36,H 28, 39 



Phase transitions [3 
equation similar to (J2J): 



40. 4 



42 



43 



44 1 can be described by an 



4>tt + 70* 



dU{4>) 



+ F(x,y,z,t). 



(15) 



Recall that {7(0) possesses two minima (0i and ^3) and a maximum (0! < 2 < 03)- 

Suppose for the moment that f7(0i) < U(<ps). The main phenomenon for F = is that 



of "nucleation" where drops or bubb 



33, 



34 



35 



36 



37 



38, 



39 



40 



41 



42 



es grow or disappear depending on their Gibbs energies 



43, 



44). 



The relevant issue is the transition from the metastable state 03 to the state 0i. In phase- 
transition theory, researchers consider the existence of a critical "germ" for the transition. 
Field configurations 0(x, t) with a "radius" larger than that of the critical germ should grow 
in order to produce the transition to the state <\>\. 

Let us first discuss this situation in the context of the one-dimensional equation 

dU(<f>) 



(16) 



For i7(0i ) < {7(03), it is possible to show (see Refs. |33j and |45|) that there exists an 
unstable stationary critical solution that plays the role of a critical germ. This solution is 
represented in Fig. H21 
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FIG. 12: A "bubble" formed by a kink-antikink pair in Eq. (|T7|) 



This solution can be interpreted as a kink-antikink pair 49, 46]. 

In general, a kink and an antikink attract each other 4(1 L^. On the other hand, a 
constant external force acts on a kink and an antikink in different directions. So the solution 
represented in Fig.^Jis possible when the kink and the antikink are at some critical distance 



where the attraction force and the external force balance each other 
The critical solution is unstable 



3, 



46, 4 



461 ] . Perturbations of the critical solution can lead to 
the growth of the critical "bubble" initiating the transition to the state (j) — "Bubbles" 
with a radius smaller than the radius of the critical solution will disappear, whereas "bub- 
bles" with a radius larger than the radius of the critical solution will grow until the whole 
system is in the state = 0i. 

An example can be studied in the equation 

1 . 1 



4>u + i<Pt - 02 



-0+ -( 
T 2 



F , 



(17) 



where F = const. 

Note that Eq. (pH can be written in the form (|T6|) . where 

f/(0) = i(0 2 -l) 2 -F o 0. 



(18) 



If F < 0, then U{4>i) < U((f) 3 ). Eq. (|T7|) possesses a bell-shaped stationary solution as 
that shown in Fig. El This solution is unstable j^Q). Similar but different bell-shaped 
initial states will evolve either to the state = 0i or the state = 03. If the kink and 
the antikink that form the one-dimensional "bubble" are too close, then the bubble will 
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FIG. 13: Inhomogeneous external perturbation able to stabilize the bubble in Eq. (|19|) 



collapse. If the distance between the kink and the antikink is larger than the "width" of the 
critical germ, then the "bubble" will grow. 
If F = 0, all the initial bubbles will collapse. 

Can we stabilize the bubbles of phase = 0i inside the "sea" of pase = 03 using 
external perturbations? 

Let us analyze what happens when a kink-antikink pair as that shown in Fig. ^] is in 
the presence of an external perturbation like that shown in Fig. EH 

Note that these graphs look similar but they represent very different physical quantities. 

If the kink of Fig. EDis close to the "right" zero of function F(x), it will "feel" an effective 
potential well (see Fig. I14jl. The same will happen to the antikink if it is close to the "left" 
zero of F(x) in Fig. H21 (See the potential in Fig. (|I5j).) 

If the potential barrier is sufficiently high, the attraction force between the kink and 
the antikink will not lead to the "bubble" collapse. The attraction force decays rapidly 
(exponentially) with the distance, so for large "bubbles" practically "any" F(x) with the 
shape of Fig. [13] is good for stabilizing the "bubble" . 

On the other hand, even small bubbles can be stabilized using an external perturbation 



X = Xl,2 



F(x) with two zeroes x\ and x 2 (as in Fig. IT^j) such that f ^jjjp 
That is, the change of sign in F(x) should be sufficiently sharp. 

We have checked numerically the existence of the stable "bubbles" 

Let us experiment with the following equation 

1 !,3 



is sufficiently large. 



0« + l<t>t - 0,x - 2^ + 2^ = ( 19 ) 
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FIG. 14: Effective potential for the kink that forms the "bubble" shown in Fig. (|12|) in the presence 
of perturbation shown in Fig. (jlMjl 



FIG. 15: Effective potential for the antikink that forms the "bubbles" shown in Fig. (|12|) in the 
presence of perturbation shown in Fig. (|13|) 

where 

F(x) = A {tanh [B(x - d)} - tanh [B(x + d)} + e} (20) 

For instance, we have established that the perturbations (}2Tjj) . with the parameters A = 1, 
B = 0.2, d = 5, e = 0.7, can sustain a stable bubble. 

An initial configuration for <p(x, t) where the kink and the antikink are not exactly in 
the zeroes of F(x) will, anyway, tend to the asymptotically stable bubble. That is, there is 
a bubble which is a spatiotemporal attractor. All close initial configurations will tend (for 
t — > oo), to the asymptotically stable bubble. 

This phenomenon occurs in the case that a bubble is considered as an initial configuration. 
That is, a bubble had been created before the "force" F(x) is applied. However, it is 

16 



interesting to remark that for sufficiently large A and B, the perturbation itself can create 
the bubble. That is, there is no need for a preexisting bubble. The bubble is created by the 
perturbations and, then, the perturbation can sustain the stabilized bubble. 

If A is too small, it is possible that F(x) cannot sustain a bubble. An example is produced 
with the following parameters A = 0.25, B = 1.0, d = 1.0, e = 0.8. All initial conditions 
lead to the bubble collapse. 

Once we have an external perturbation F(x) able to sustain a stable bubble, we can 
transport the bubble if we can move F(x) as a whole. This phenomenon is in the same 
spirit of all the theory discussed above concerning the soliton transport. 

In the case of the bubble, as in many of the discussed situations with solitons, there is a 
limit for the velocity of the perturbation as a whole. When the shape represented in Fig. 
is moved too fast, the bubble can be left behind, and then it will collapse. 



V. THREE DIMENSIONAL BUBBLES 



What about the bubbles in three dimensions? 

Let us concentrate on solutions with radial symmetry. 

The stationary solutions of Eq. (|T5|) (F(x,y, z,t) = 0), satisfy the equation 

Ar + -A- d -^ = 0. (21) 

r d(p 

It is possible to prove that (for U(4>i) < U (03)) there exists a solution with a minimum 
in r = and that asymptotically it tends to 03 as r — > oo. This solution corresponds to the 
critical germ discussed above. Again as before, this critical solution is unstable. 

In general, exact solutions for the three-dimensional case of the nonlinear equation (|21|) 
are very difficult to find . Nevertheless, we have developed an inverse-problem approach 
that allows us to find exact solutions to a class of topologically equivalent systems. 

For instance, a wealth of work [33] has been dedicated to find a solution to equations of 
type 

rr + -0 r - B 2 <p + D 2 3 + T(0) = 0, (22) 

where T(0) possesses only terms of higher order than 3. This equation is important in field 
theory. 
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We have found that function 

Atcmh(Br) 



(23) 
r cosh(5r) 

is the exact solution to equation rr + 2 r — B 2 (p + ^p0 3 + T((f>) = 0, where T(0) contains 
only terms of order higher than 3. 

Likewise, it can be proved that function 

Ataiah(Br) , . 

r cosh(i3r + D) 

is the critical germ solution of an equation of type (|21|). where £7(0) possesses all the prop- 
erties that we need: two minima and a maximum, U(<f)i) < £7(03), etc. 

As in the one-dimensional case, small bubbles will collapse, and bubbles with a radius 
larger than that of the critical germ will grow. 

In the context of the three-dimensional equation: 

0« + 70*-V 2 0-^+^0 3 = O, (25) 

all the bubbles are non-stationary. However, we should say that very large bubbles are 
quasi-stationary because the interaction between the walls is so small that (although they 
should collapse eventually) their life-time is very large. 
What about inhomogeneous perturbations as in Eq. (jl5|) ? 

Suppose we have an initial configuration 0(x, t = 0) with the shape of a spherical bubble 
of phase — 0i inside a "sea" of phase 03. This means that there is a spherical spatial 
sector where ~ 0i, and, outside this sector, w 3 . Space zones where F(x, y, z) < will 
"try" to expand the bubble. Space zones where F(x, y, z) > will "try" to make the bubble 
collapse. 

In fact, as in the one-dimensional case, we can stabilize a spherical bubble using an 
inhomogeneous external perturbation F(x,y,z) such that F(x,y,z) < for points (x,y,z) 
that are inside the bubble and F(x, y, z) > for points (x, y, z) outside the bubble. 

To ensure the stability of the bubble, the sign change of F(x, y, z) should be sharp. We 
have an exact solution. 

Let us define the equation 

<f>tt + l<f>t - V 2 - ^0 + ^0 3 = F(r), (26) 
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where 



F(r) 



2 

— tanh 



A (A 2 - 1) | tanh 



A 



r-R) 



+ 2 { (0i + <h) [0i02 + A (0i + 2 ) + A 2 ] } 



lr + 02r 



-A+-A\ 
2 2 ' 



r = ^x 1 + y 2 + z 2 , 0i = A tanh [£(r - i?)], 2 = -A tanh [S(r + R)}; A, B and are real 
parameters. 

Here F(r) has the desired property that F(r) < in a spherical sector around r = 0, and 
F(r) > outside this spherical sector. 
There is an exact bubble solution to Eq. 



A < tanh 



^(r-R) 



tanh 



^(r + R) 



(27) 



But we can say more. For R > 5, A > 1, this bubble solution is stable. Manipulating 
F(r) we can control this bubble. In fact, we can move this bubble to another point in space. 



VI. BUBBLES IN PHASE TRANSITIONS 



Many works have been dedicated to the nucleation process in phase transitions 33|, |34 , 



41 



42,0, 44,Q 



Here we wish to focus on the works based on the 
42, 3,Q- 



0, 0, H H 

theory of relaxation of metastable states 

In Hefsfl Q Q, the nuclear in the notable phaae nea r the c riti ca, pen, of 
a thermodynamic system is described as a relaxation of a metastable state of the order- 
parameter field. 

A metastable state of matter can be produced in a first-order phase transition. Its 
relaxation into a thermodynamically stable state is facilitated by the formation of a critical 
nucleus of a new phase. 

The system is described by a scalar field 0(x, t) for which a relaxation equation is obtained: 



1 dcj) _8H 

r n m ~~50 +text1 



(28) 
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where H[<f>] is the system free-energy which is taken in the Landau form: 

H\<P] = \j {c ( V0) 2 + ^ 2 + l -g<p A - 2/10} dx. (29) 

Here T n is a kinetic coefficient and f ext is an external perturbation. 

At h = and /J = \i c < 0, there is a critical point in the system. The line h = (/x < /i c ) 
is a line of first-order phase transition. 

The Eq. f!28j) is an overdamped version of the 4 equation discussed above. This equation 
can be written in dimensionless units as our equations (j2J). The solution we have been 
investigating describes a nucleus of phase (0) = —1 in phase (<p) = 1. At h < 0, the nuclei 
with small radius attenuate, while the nuclei with large radius increase. The phase ((f)) = 1 
at h < is metastable-unstable to formation of a nucleus with a large radius. 



The parameters of the system ()28H29jl for the liquid-vapor first-order phase transition 
depend on the thermodynamics quantities, e.g., pressure and temperature 3,13- 

If an external perturbation (say a laser beam) creates the thermodynamics conditions for 
a first-order transition (e.g., h < h c < in the context of system (|28H29| ) in a small zone 
of a liquid mass, then a vapor bubble could be generated and sustained. 

We should say here that in the framework of equation 

j<f> t - V 2 - 1 -<P + ^ = F(x, y, z), (30) 

in order to create a bubble able to increase to a macroscopic radius, F(x, y, z) must take 
negative values in a limited spatial zone in such a way that F(x,y,z) < — On the 
other hand, once the bubble has been generated, the perturbation F(x, y, z) necessary to 
sustain the bubble should be "negative" inside the bubble and positive "outside", but the 
absolute value of F(x, y, z) in the "negative" part should not be so high. In other words, the 
"intensity" of the perturbation needed for generating the bubble is larger than the needed 
"intensity" for stabilizing the bubble. 

VII. REAL EXPERIMENTS: BUBBLES INDUCED AND CONTROLLED BY 
LASERS 

Now we will discuss several experiments about bubble formation in a weak absorbing 
liquid and its subsequent trapping under the action of a cw laser beam of relatively low 
power. 
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FIG. 16: Experimental Setup 



The generation of bubbles by intense light pulses was studied in the papers [48, 49]. On the 
other hand, the pioneering works about laser bubble trapping can be found in Refs. I50L [51 1 . 
Here we present a review of previous 48, |^ and very recent experiments concerning new 



effects related to laser bubble trapping, in the light of the theory we are developing. 

High-intensity light induces a local first-order liquid-vapor phase transition. If light 
power is reduced after the bubble formation, the bubble remains trapped. If the intensity 
is very low then the bubble collapses. 

The experimental setup is shown in Fig. EH 

In Figure El w e show the experimental set-up used for the laser bubble generation and 
subsequent trapping. A CW 200 mW diode pumped green laser (1=532 nm) is directed 
vertically down toward a horizontally located sample cell using the mirror Ml. A 15-cm 
focal length lens L focuses this light onto the sample. By moving this lens we can change 
the radius of the focal spot. The sample is a 4- mm thick glass cell containing an iodine 
ethanol solution of high concentration (1 to 10 ml). The focused beam heats this solution, 
generating locally a small bubble at power levels larger than 60 mW. Following the light 
absorption a distribution of temperatures is induced in the focal region and the bubble 
surface. Because the surface tension depends upon the temperature, a spatial distribution 
of surface tension is generated on the bubble. The generated gradient of surface tension 
induces a force that traps the bubble at the point of maximal temperature (focal point). 
This generated force together with the floatation Archimedes force make the bubble sticks 
to the upper surface of the cell. By moving the laser beam using the mirror Ml, the bubble 
follows the beam spot demonstrating the actual bubble trapping by the laser beam. The 
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FIG. 17: Trapped bubble using a laser. The laser power is 1.7 mW (left), 0.5 mW (center) and 
0.15 mW (right). The bubble's diameters depend on the laser power. 

bubble can be observed using the mirror M2 and a CCD camera connected to a computer. 

The used liquids were ethanol, benzene, and ethylene glycol. The liquids were colored 
with iodine or a dye in order to increase absorption. 

We will focus here on the bubble behavior. There are several interference phenomena 
that will not be discussed here. 

In Fig. llTHT^l pictures of the trapped bubbles are shown. 

In one of the experiments the bubble is generated at a laser power of 70 mW. The bubble 
size decreases with decreasing power. 

These phenomena can be described in terms of thermal gradients and forces 5l|. However, 



we wish to stress the universal character of the phenomena related to bubble formation and 
transport. These effects can occur in many other physical systems with phase transitions, 
solitons, domain walls and bubbles of one phase inside other phase, 
The bubble can be controlled by the laser. 

If after creating the bubble, the light beam is defocused in such a way that its spot 
size becomes much larger than the radius of the generated bubble, then changes of the 
bubble radius occur until this radius is again stabilized. This can be explained using the 
interpretation discussed above where the bubble is described as trapped inside a stable 
equilibrium position or potential well. 

The bubbles are created with a laser power of 70 mW. Then the power is decreased about 
10 times. 

In Fig. El we can observe a trapped bubble with a diameter 0.4 mm using a laser power 
1.7 mW. When the laser power is 0.15 mW the bubble diameter is 0.2 mm. 

Fig. fl8l shows the transport of the bubble using the laser beam. With the help of a mirror 
the light beam is moved. The trapped bubble can be moved as the hight beam is moved 
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FIG. 18: Transport of the bubble: (left) The bubble is shown at its initial position, (right) The 
bubble is moved to the left superior corner 











• 

















FIG. 19: The bubble can collapse when the laser power is small. 



slowly. If the motion of the light beam is too fast, then the bubble is left behind and then 
it collapses. 

In Fig. (left), the bubble is shown at is initial position. Then it is moved to the left 
superior corner using the laser beam. 

Fig. El shows the collapse of the bubble while the light power is maintained in 0.1 mW. 



VIII. CONCLUSIONS 



We have investigated generalized nonlinear Klein-Gordon equations perturbed by moving 
inhomogeneous external perturbations. 

We have shown that a kink-soliton can be transported by an external wave, provided the 
shape of the wave and its velocity satisfy certain conditions. For instance, a moving external 

perturbation F(x — vt), where function F(y) possesses a zero y such that |p > 0, can 

v y=yo 

carry a kink-soliton provided its velocity does not exceed certain limit. Figures ^ H3 and 
El show examples of these transporting inhomogeneous perturbations. Sinusoidal waves as 
that shown in Fig. 03 can also be successful in transporting a kink-soliton under certain 
conditions. 

We have also determined the necessary conditions to be satisfied by an external pertur- 
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bation to stabilize a bubble of one phase inside another. 

If this stabilizing inhomogeneous external perturbation is moved, then the bubble can be 
transported to another space point. 

We have presented the results of real experiments with laser-induced vapor bubbles in 
liquids. The bubbles can be created, trapped, stabilized and transported by a laser beam of 
relatively low power. 

The investigated equations are very general. Besides, the concepts and results that lead 
to the existence of the reported phenomena are also universal. So we believe that these 
phenomena can be observed in many other physical systems. For instance, the stabilization 
and transport of bubbles using an external perturbation can occur in other systems under- 
going a phase-transition. So the mentioned bubbles (or drops) can be superconducting, 
ferromagnetic, superfluid, etc. 

On the other hand, the discussed transport of solitons and kink-antikink bubbles can 
be very important in spin systems, molecular motors, charge-density waves, ferroelectric 
materials and Josephson junctions. In particular, in long Josephson junctions the kink- 
solitons are ffuxons. These fluxons can be controlled using external currents [3, [yj which 
can be very relevant in new communication and computer technologies where the fluxons 
are used as very stable information units. 

Moreover, the kinks are examples of topological defects. Many results about the kinks 
can be generalized to other topological defects as vortices and spiral waves which are being 
studied intensively in countless systems nowadays [3,0, 
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